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ABSTRACT 

We examine the origin and evolution of correlations between properties of supermassive black holes (BHs) 
and their host galaxies using simulations of major galaxy mergers, including the effects of gas dissipation, cool- 
ing, star formation, and BH accretion and feedback. We demonstrate that the simulations predict the existence 
of aBH "fundamental plane" (BHFP), of the formM BH oc o-3-Oio.3flO.43io.i9 orM BH ocM°- 54±0 - 17 o- 2 - 2±a5 , sim- 
ilar to relations found observationally. The simulations indicate that the BHFP can be understood roughly as a 
tilted intrinsic correlation between BH mass and spheroid binding energy, or the condition for feedback cou- 
pling to power a pressure-driven outflow. While changes in halo circular velocity, merger orbital parameters, 
progenitor disk redshifts and gas fractions, ISM gas pressurization, and other parameters can drive changes in 
e.g. a at fixed M*, and therefore changes in the M Bll -a or M BH -M* relations, the BHFP is robust. Given 
the empirical trend of decreasing R e for a given M* at high redshift (i.e. increasingly deep potential wells), the 
BHFP predicts that BHs will be more massive at fixed M*, in good agreement with recent observations. This 
evolution in the structural properties of merger remnants, to smaller R e and larger a (and therefore larger Mbh, 
conserving the BHFP) at a given M*, is driven by the fact that disks (merger progenitors) have characteristically 
larger gas fractions at high redshifts. Adopting the observed evolution of disk gas fractions with redshift, our 
simulations predict the observed trends in both R e (M*) and Mbh(M»). The existence of this BHFP also has 
important implications for the masses of the very largest black holes, and immediately resolves several appar- 
ent conflicts between the BH masses expected and measured for outliers in both the M^-a and Mbh~M* 
relations. 

Subject headings: quasars: general — galaxies: active — galaxies: evolution — cosmology: theory 



1. INTRODUCTION 

Correlations between the masses of supermassive black 
holes (BHs) in the centers of galaxies and the properties 
of their host spheroids (e.g., Kor mendy & Richstonelll995h 
imply a fundamental bond between the growth of BHs 
and galaxy formation. A variety of such correlations 
have been identified, Unking B H mass to host luminosity 
dKormendv & Richstond 119951). mass (|Magorrian et alj 
19981). veloci t y dis persion dFerrarese & MerrittI l2000t 
Gebhardt et alj 12000]). concentration or Sersic index 
dGraham et al J 120011: iGraham & Driver! 120061) . and binding 
energy dAller & Richstone1l2007l) . However, the connection 
between these relationships is obscured by the fact that the 
properti es of the host sp heroids are themselves correlated 
(see e.g. iNovak et al.ll2.006l for a comparison of the observed 
relations). The lack of a clear motivation for favoring one 
correlation over another has led to considerable debate over 
the interpretation of systems deviating from the mean corre- 
lation between host properti es, and over the demographics of 
the most massive BHs (e.gjBernardi et al. 2006; Lau er et all 
2006b; iBatcheldor et al. 2006: IWvithell2006l). 

Analytical estima tes (e.g., ISilk & Reesl 119981; 
iBurkert & Silkl 1200 lh and studi e s usin g simulations 
dCox et al l 1200 6b; Rob ertson et al.l l2006bl) have demon- 
strated that these corr elations, in particul ar the Mbh ~ Q" 
dPi Matteo et al.l 120051) and M BH - M* dRobertson et al.l 
2006c) relations, can be reproduced in feedback-regulated 
models of BH growth. However, determining the fundamen- 
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tal character and evolution of these correlation s with r edshift 
is critical for informi ng analytical models ( e.g., Croton l2006l) 
and simulations jP i Matteo et al.l 120051; Robertso n et al.l 
l2006d:lHopkins et alj2005lf which follow the co-formation of 
BHs and bulges, as well as theories which relate the evolution 
and statistics of BH formatio n and q uasar ac tivity to galaxy 
mergers (e.g., iHopkins et all l2006al l2007fl l2006d l2007cf) 
and to the r emnant spheroid population (Hopki ns et al.l 
l2006cl l2007bl) . Likewise, the significance of observations 
tracin g the buildup of spheroid populations (e.g JCowie et al.l 
1996) and associations between spheroids in formation, 
mergers, and quasar hosts (Hopkins et al. 2007a) depends on 
understanding the evolution of BH/host correlations. 

Unfortunately, efforts to directly infer these correlations at 
redshifts z ^> are difficult and still limited by the small num- 
ber of observable hosts. Furthermore, without understand- 
ing the fundamental nature of the correlations between BH 
and host properties, it is difficult to interpret these observa- 
tions, as they do not all probe the same correlations. Conse- 
quently, different groups have arrived at seemingly contradic- 
tory conclusions. Velocit y dispersion measu rements have fa- 
vored both no evolution ( Shields et al. 2003, from QUI veloc- 
ity dispersions) and substantial evolution dShields et alJl2006t 
IWoo et aii r2006, from CO dispersions and spectral template 
fitting') . BH clustering measurements (jA delberger & Steidel 
l200l lWvithe & Loeb 2005; Hopk ins et alJl2007et iLidz eTal 



2006) suggest moderate evolution in the ratio of BH to host 
halo mass at redshif ts z ~ 1 - 3. Pir ect host 7?-band luminos- 
ity measurements dPeng et al.ll2006l) and indire ct comparison 
of qu asar luminosity and stellar mass densities ( Merloni et al.1 
120041) or BH and stellar mass functions (Hopkins et al. 2006d) 
similarly favor moderate evolution in the ratio of BH to 
host spheroid stellar mass occurring at z > 1, and dynami- 
cal masses from CO measurements suggest that this evolution 
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may extend to z ~ 6 dWalter et alj |2004). Better understand- 
ing the dependence of BH mass on host properties can provide 
both a self-consistent paradigm in which to interpret these ob- 
servations and (potentially) a physically and observationally 
motivated prediction of their evolution. 

One possibility is that these different correlations are pro- 
jections of the same "fundamental plane" (FP) relating BH 
mass with two or more spheroid properties such as stellar 
mass, velocity dispersion, or effective radius, in analogy to 
the well-established fundamental plane of spheroids. For the 
case of spheroids, it is now understood that various correla- 
tions, including the Faber-Jackson relation (Faber & J ackson! 
119761) between luminosit y (or effectively ste llar mass M») and 
velocity dispersion a, the Kormendy] (1 19771) relation between 
effective radius R e and surface bright ness I e , and the s ize- 
luminosity or size-mass relations (e.g., She n et aLll2003l) be- 
tween R e and M„, a re all projections of a fundamental plane 
relatin g R e oc a a /f dDressleret al.ll 19871: iDiorgovski & Davis! 
f!987h . 

In their analysis of the relation bet ween BH mass and hos t 
luminosity or dynamical mass, Md v n , iMarconi & Hunt! (120031) 
(see also Ide Francesco et al.l 12006b noted th at the residuals 
of the Mbh-ct relation (effectively Mbh/c 4 ; Tremai ne et al.l 
l2002h were significantly correlated with the effective radii 
of the system s in thei r sample. A more detailed study by 
Hopkins et al. (2007d) (henceforth Paper I) found that the 
observations indeed favor a fundamental plane-type relation 
between Mbh and any two of R e , a, or M*, at > 3<r con- 
fidence. This observed, low redshift BHFP has the form 
M BH oc CT 3.o±o.3^o.43±o.i9 or Mm K M o.54±o.iv CT 2.2±o. 5> with 

the early-type FP providing a tight mapping between M*, R e , 
and a. Given the mean correlations between e.g. R e and M* or 
M, and a, the previously recognized correlations between BH 
mass and spheroid mass, luminosity, velocity dispersion, and 
binding energy can all naturally be explained as projections 
of this intrinsic BHFP. However, these observations remain 
limited in the range of systems they probe, and are restricted 
to passive, local spheroids. We therefore wish both to under- 
stand the origin of this BHFP and how it should (or should 
not) evolve with redshift. 

In this paper, we investigate the nature of the correlation 
between BH mass and host properties and the existence of a 
fundamental plane relating BH mass and spheroid mass, ve- 
locity dispersion, and effective radius. In § [3] we describe a 
large set of numerical simulations which we use to study and 
predict the nature of the BH-host correlations under a wide 
variety of conditions, and in §|2]we describe the observational 
data sets we compile to study the observed correlations. In 
§ |4] we describe the correlations determined from both sim- 
ulations and observations, and then analyze the correlations 
between residuals in e.g. the Mbh-ct relation and seco ndary 
properties such as R e and M*, which leads us in § !4.2l to dis- 
cuss the fundamental plane relating BH mass and a, R e , and 
M». § |5] discusses the implications of this relation for pre- 
dicting BH masses and demographics, and § [6] considers the 
physical origin of the BHFP relation. In § [7] we study how 
various theoretical quantities or initial conditions drive sys- 
tems along the BHFP relation and, as a consequence, drive 
evolution in the various projections of the BHFP, and in § [8] 
apply this to understand the observed evolution with redshift 
in the Mbh~M» and Mbh-ct relations. We summarize our 
conclusions and discuss future tests of our proposed relations 
in §|9] 

Throughout, we adopt a J7m = 0.3, ft\ = 0.7, Hq = 



70km s" 1 Mpc" 1 cosmology (and correct all observations ac- 
cordingly), but note this choice has little effect on our conclu- 
sions. 

2. THE DATA 

The observational data set with which we compare is de- 
scribed in detail in Paper I, but we summarize it here. We 
consider the sample of local BHs for which masses have been 
reliably determined via either kinematic or maser measure- 
ments. Specifically, we adopt the sample of 38 local systems 
for which values of Mbh, o, Re, M d yn , and bul ge luminosities 
are co mpiled in Marconi & Hunt! (120031) and iHE ring & Rix 
(120041) (see also iMagorrian et al.l 11998b iMerritt & Ferraresel 
l2001tlTremaine et al.ll2002l) . We adopt the dyn amical masses 
from the more detailed Jeans modeling in Haring "& Rixl 
(2004). We estimate the total stellar mass M* f r om th e to- 
tal K-b&nd luminosity given in IMarconi & Hunt! (120031) . us- 
ing the Zr- band mass-to-ligh t ratios as a function of lumi- 
nosity from lBell et al.l d2003l) (specifically assuming a "diet" 
Salpeter IMF, although this only affects the absolute normal- 
ization of the relevant relations). In Paper I it is also noted 
that the choice of these mass-to-light ratios as opposed to 
those determined from e.g. photometric fitting makes little 
diffe rence. We take measure ments of the Sersic index n s 
from [Graham &~Driverl 12006). Where possible, we update 
measurements of R e , a and n s with more recent v alues from 
lLauer et al.l (120051 l2006al) : iMcDermid eTail (120061) and from 
iKormendy et al.l(l2007l) . which extends the baseline of surface 
brightness measurements allowing more robust estimates of 
n s and/? e . 

Although it should only affect the normalization of the rela- 
tions herein, we note that our adopted cosmology is identical 
to that used to determine all quoted values in these works. The 
concentration index R30/R50 for the observed systems is cal- 
culated assuming a Sersic profile with the best-fit n s . When 
we fit the observations to e.g. the mean Mbh-c relation and 
other BH- host relations, we consid er only the subsample of 27 
objects in Marc oni & Hunt! (120031) which are deemed to have 
'secure' BH and bulge measurements (i.e. for which the BH 
sphere of influence is clearly resolved, the bulge profile can 
be well-measured, and maser spots (where used to measure 
Mbh) are in Keplerian orbits). Our results are not qualita- 
tively changed if we consider the entire sample in these fits, 
but their statistical significance is somewhat reduced. 

3. SIMULATIONS 

3.1. Methodology 

Our simulations are taken from Ro bertson et al.l (|2P06b), 
who utilize a set of several hundred simulations of major 
mergers to study the properties of remnants on the early-type 
galaxy FP. The properties of the models are discussed in de- 
tail therein, but we briefly review them here. The simulations 
were performed with the parallel TreeSPH code GADGET-2 
(Springel 2005) , based on a fully conservative formulation 
( Springel & Hernquist 2002) of smoothed particle hydrody- 
namics (SPH), which conserves energy and entropy simul- 
taneously even when smoothing lengths ev olve adaptively 
(see, e.g.. lHernquislll993t IcFShea et al J 120051) . Our simula- 
tions account f or radiative cool i ng, heating by a UV back- 
ground (as in iKatz et al.l 1 19961; iDave et al.l 1 19991) . and in- 
corporate a sub-resolution model of a multiphase interstel- 
lar mediu m (ISM) to describe st ar for mation and sup ernova 
feedback (ISpringel & Hernquisll l200l ISpringel et ail 120051) . 
Feedback from supernovae is captured in this sub-resolution 
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model through an effective equation of state for star-forming 
gas, enabling u s to stably evolve disks w i th arbitrary gas 
fractio ns (see [Springel & Hernquisj] 120051; iRobertson et al.l 
l2006al) . This feedback prescription can be adjusted between 
an isothermal gas with effective temperature of 10 4 K and our 
full multiphase model with an effective temperature ~ 10 5 K. 

Supermassive black holes are represented by "sink" par- 
ticles that accrete gas at a rate M estimated from the local 
gas density and sound speed using an Eddington-limited pre- 
scription based on Bondi-Hoyle-Lyttleton accretion theory. 
The bolometric luminosity of the black hole is taken to be 
Lboi = e,Mc 2 , where e r = 0.1 is the radiative efficiency. We 
assume that a small fraction (typically w 5%) of Lboi couples 
dynamically to the surrounding gas, and that this feedback is 
injected into the gas as thermal energy, weighted by the SPH 
smoothing kernel. This fraction is a free parameter: we adjust 
it to match the normalization o f the local Msh ~ o r elation as 
in lDi Matteo eTaII(l2005ll2007h : ISijacki et~ai1(l2007l) . We em- 
phasize that this controls only the normalization of this rela- 
tion; i.e. inefficient feedback coupling means a BH must grow 
proportionally larger in order to couple the same energy to the 
ISM and self-regulate, but the scalings of BH mass with a and 
host properties (i.e. slopes of the BH-host relations and corre- 
lations between residuals in these relations) are not changed. 
Because our comparisons throughout are based on the relative 
scalings of BH mass with host properties, this normalization 
choice is simply a matter of convenience. For now, we do not 
resolve the small-scale dynamics of the gas in the immediate 
vicinity of the black hole, but assume that the time-averaged 
accretion rate can be estimated from the gas properties on the 
scale of our spatial resolution (roughly w 20 pc, in the best 
cases). 

The progenitor galaxies in the mergers are constructed fol- 
lowing Springe l et al.l (120051) . For each simulation, we gen- 
erate two stable, isolated disk g alaxies, ea c h with an ex- 
tended dark matter halo with a iHernquisll d 1990b profile, 
motivated by cosmolog ical simulations (e.g.. iNavarro et al] 
ll996tlBusha et al.l2005l) . an exponential disk of gas and stars, 
and (optionally) a bulge. The galaxies have total masses 
M v i r = V v 3 ir /(10G//o) for z = 0, with the baryonic disk having 
a mass fraction m<j = 0.041, the bulge (when present) hav- 
ing nvo = 0.0136, and the rest of the mass in dark matter. 
The dark matte r halos are assigned a c oncentration parame- 
ter scaled as in IRobertson et al. l d2006d) appropriately f or the 
galaxy mass and redshift following Bullock et al. (2001). The 
disk scale-length is computed based on an assumed spin pa- 
rameter A = 0.033, chosen to be n ear the mode in the A dis- 
tribution measured in simulations (IVitvitska et alJl2002h . and 
the scale-length of the bulge is set to 0.2 times this. 

Typically, each galaxy initially consists of 600,000 dark 
matter halo particles, 20,000 bulge particles (when present), 
40,000 gas and 40,000 stellar disk particles, and one BH pani- 
cle. We vary the numerical resolution, with many simulations 
using twice, and a subset up to 128 times, as many particles. 
We choose the initial seed mass of the black hole either in 
accord with the observed Mbh-c relation or to be sufficiently 
small that its presence will not have an immediate dynamical 
effect, but we have varied the seed mass to identify any sys- 
tematic dependencies. Given the particle numbers employed, 
the dark matter, gas, and star particles are all of roughly equal 
mass, and central cusps in the dark mat ter and bulge are re a- 
sonably well resolved (see Figure 2 in Spri ngel et al.l l2005). 

We consider the set of several hundred simulations from 



IRobertson et al.l d2006bl) . in which we vary the numerical res- 
olution, the orbit of the encounter (disk inclinations, peri- 
center separation), the masses and structural properties of 
the merging galaxies, initial gas fractions, halo concentra- 
tions, the parameters describing star formation and feed- 
back from supernovae and black hole growth, and initial 
black hole masses. The detailed list of v a ried pr operties is 
given in Tables 1 & 2 of IRobertson et al] ([2006b). For ex- 
ample, the progenitor galaxies have virial velocities V v j r = 
50,80, 115, 160,226,320, and 500 kms" 1 , and are constructed 
to match disks at redshifts z = 0,2, 3, and 6, and our simu- 
lations span a range in final black hole mass Mbh ~ 10 5 - 
1O 1O M0. The extensive range of conditions probed pro- 
vides a large dynamic range, with final spheroid masses 
spanning M* ~ 10 8 - 1O 13 M , covering the entire range of 
the observations we consider at all redshifts, and allows 
us to identify any systematic dependencies in our models. 
We consider initial disk gas fractions (by mass) of / gas = 
0.05, 0.2, 0.4, 0.6, 0.8, and 1.0 for several choices of virial 
velocities, redshifts, and ISM equations of state. 

The results described in this paper are based primarily on 
simulations of equal-mass mergers; however, by examining 
a small set of simulations of unequal mass mergers, we find 
that the behavior does not change significantly for mass ratios 
down to about 1:3 or 1:4, below which mass ratio the mergers 
produce neither substantial BH nor bulge growth, and there- 
fore the are no longer appropriate to compare to local relations 
between BHs and massive spheroids. 

3.2. Analysis 

Each simulation is evolved until the merger is complete 
and the remnants are fully relaxed, typically ~ 1 — 2 Gyr 
after the final merger and coalescence of the BHs. We 
then measure kinematic properties of the re mnants following 
IRobertson etail (l2006bh : ICo"x et alJ (l2006bl) . The effective ra- 
dius R e is the projected half-mass stellar effective radius, and 
the velocity dispersion a is the average one-dimensional ve- 
locity dispersion within a circular aperture of radius R e . Pro- 
jected quantities such as R e , a, and the stellar surface mass 
density I e = M^{r < R e )/TrR 2 e are averaged over 100 random 
lines of sight to the remnant. Throughout, the stellar mass M* 
refers to the total stellar mass of the galaxy, and the dynamical 
mass Mdy n refers to the traditional dynamical mass estimator 

M dyn = k—^-, (1) 

where we adopt k = 8/3 (although this choice is irrelevant as 
long as we apply it uniformly to both observations and sim- 
ulations). We define cj> c as the gravitational potential from 
the galaxy (excluding the BH itself) at the location of the 
BH (given the typical resolution over which this quantity is 
smoothed, this is effectively the conve rged bulge pote ntial at 
r = 0, in the absence of the BH; see ISpringell 12005). As a 
concentration index we adopt the ratio of half-mass radius 
R e = R$o to 30%-mass radius R30, and measure n s from pro- 
jected mock images following Krause et al. (2007, in prepa- 
ration). We note that, for convenience, / gas typically refers to 
the gas fraction in the merging disks when the simulations are 
initialized, but show in §|7]that our results are unchanged (al- 
though / gas itself systematically shifts) regardless of the time 
before the merger at which we choose to define the gas frac- 
tion of the systems. 

4. THE LOCAL BH-HOST CORRELATIONS 
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a [ km s' 1 ] log{ M, / M e ) 

FIG. 1. — Location of our simulat ion merger remnant spheroids on the 
Mbh~ c and Mbh — M* relations (as in Di Matteo et al. 2005 ; Robertson et al. 
2006c). Solid lines show the observed relations from lTremaine et al.l 420020 
and Haring & Rix ( 2004). Error bars in a show the dispersion across 100 ran- 
dom viewing angles. The simulations reproduce well the observed relations 
over a wide dynamic range. As discussed in the text, there are a number of 
likely reasons for the slight (0.2 dex) normalization offset in M» , but insofar 
as the slopes of the relations are identical, this has no effect on our analysis. 

4.1. One-to-One Relationships 

Figure Q] shows the location of our simulation remnants on 
the the Mbh ~ <J and M bh ~ M, relations. As demonstrated by 
iDi Matteo et al.l (120051) . they agree well with the observed re- 
lations over a large dynamic range. Critically, although mod- 
ifying our feedback prescriptions and, as we show below, ad- 
justing the kinematic properties of the remnants by changing 
e.g. orbital parameters and gas fractions of the merging sys- 
tems can shift the normalization of the relations, the slopes 
are not adjustable or tunable, but are a natural consequence of 
self-regulated BH growth. We note that there is a slight off- 
set between the normalization (but not sl ope) of our predicte d 
M B h - Af* relation and that measured by Hari ng & Rixl d2004l) . 
but a weaker offset in Mbh _ c. This owes to the fact that, at 
fixed er, our simulated systems typically have slightly larger 
(mean offset sa 0.16dex) stellar masses than the observed sys- 
tems on the Mbh _ cr relation. There are a number of possible 
explanations for this offs et in the Faber- J ackso n (M*[er]) re- 
lation. As pointed out by Bernardi et al.l d200 6). the systems 
with measured BH masses in lHaring & Rixl d2004) actually lie 
above the Faber-Jackson relation observed for typical early- 
type galaxies, perhaps owing to a selection bias. Their esti- 
mate of the magnitude of this bias is quite similar to the offset 
here and can completely account for our results. 

Furthermore, we show below that at fixed M», changing 
the gas fractions of the merging systems, their orbits, or their 
structural properties can systematically drive changes in a. 
This means that the precise normalization of the observed 
Faber-Jackson relation depends, in detail, on the exact star 
formation and merger histories of the systems observed. Since 
we do not model these cosmological histories, but rather iso- 
late different mergers in order to study how these changes are 
driven, it is not surprising that the normalizations of the re- 
lations do not happen to match perfectly. In any case, we 
are interested in how offsets or evolution from one relation 
or the other are produced, and how such residuals may scale 
with other host properties, for which the actual normalization 
of the relation factors out completely. It therefore makes no 
difference to our analysis (although we have considered both 
cases) whether we compare to our full suite or select only a 
subset of simulations which reproduce the (mean) normaliza- 
tion in the observed Faber-Jackson relation. 

Figure|2]shows the correlation between BH mass and a wide 
variety of host properties, from both our simulations and the 



observed sample. The slopes of the simulated correlations are 
essentially identical to those observed in every case. Note that 
many of the correlations are similarly tight, including the cor- 
relations with velocity dispersion a, stellar mass M», dynami- 
cal mass Mdy n , effective bulge binding energy M» a 2 , and cen- 
tral potential <p c . The best-fit correlations are listed in Table[T] 
along with the intrinsic scatter in Mbh estimated about each 
correlation from the simulations. We do not list the correla- 
tion with Mhaio, as it is clear in our simulations that BH mass 
is correlated with small-scale bulge properties (unsurprising, 
given that the central potential of the bulge is strongly baryon- 
dominated). Therefore, while there is an indirect correlation 
withMhaio through e.g. theMhaio-M, andMBH-M* relations, 
its nature depends systematically on the exact Mh a io - M* re- 
lation. 

We also note that the (relatively large) scatter in the cor- 
relations between Mbh and concentration or S ersic index in 
Figure [2] appears contrary to the conclusions of Graha met al.l 
(2001), who argue for very small in t rinsic scatter in these 
correlations. However, Nov ak et"aT1 (j2006) point out that 
uncertainty in this correlation, unlike for the Mbh - cr or 
Mbh -Md yn relations, is dominated by the measurement er- 
rors in concentration index or n s , which means that improved 
observations are needed to determine whether the relation 
is actually consistent with small intri nsic scatter. In fact , 
wh en we update the measure ments from lGraham et all (12001 1) 
and [Gr aham & Driveil (l2006l) with the n s measurements from 
iKormendv et al.l (120071) . which typically reduce the measure- 
ment error in n s from ~ 20% to < 5% (a nd in at least two 
cases 4 change n s by > 3a relative to the iGraham & Driver! 
(2006) fit), the quality of the correlation is substantially de- 
graded, and a significantly larger intrinsic scatter is implied. 

In the simulations, it is possible, given appropriate gas 
fractions or orbital parameters, to substantially change n s or 
R30/R50 at fixed stellar mass without driving a corresponding 
change in Mbh («.s also appears to be more variable sightline- 
to-sightline than other projected quantities such as a or R e ). 
As will be discussed in detail in Krause et al. (2007, in prepa- 
ration), this has the consequence that the values of n s pre- 
dicted by a uniform set of simulations are not as tightly corre- 
lated with spheroid mass as is observed (and, as a secondary 
consequence, Mbh, which is tightly correlated with M*, is 
less well-correlated with n s ). Furthermore, the weaker cor- 
relation in the simulations should not be surprising - unlike 
the M» - R e or M» - a correlations which are produced fun- 
damentally (at least to lowest order) by basic dynamical ef- 
fects (and are not especially sensitive to the exact type of 
spheroid-producing m ergers), it is commonly believed (e.g., 
IKormendv et al.l 120071) that the correlation between n s and 
M* is driven by an increasing prevalence of dissipationless 
(spheroid-spheroid) mergers at higher masses. Such merg- 
ers will, by definition, conserve the Mbh-M* relation, but 
simulations have shown that they dramatically change n s (e.g. 
iBoylan-Kolchin et al.ll2006l) and therefore the Mbh~«j rela- 
tion. If this is true, then the Mbh _ n s correlation, unlike the 
Mbh-M» relation, must fundamentally be driven by cosmo- 
logical effects (e.g. the differential contribution of different 
merger types as a function of mass and redshift), which our 
simulations do not represent. A more accurate modeling is 

4 IGraham & Driveil J2006I) quote = 3.04+°^|, 2.73^ -^, for NGC 
3377 and 4473, respectively, whereas Kormendv et al. (2007) measure n s = 
4.917, 6.04 ±0.33. The ne wer values make the mea sured Mbh discrepant 
with the Mbh _ n s relation in Graham & Drivei 1 2006) by 0.64 and 0.78 dex, 
respectively (with ss 0. 15 dex measurement errors in Mbh for each). 
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FIG. 2. — Correlations between BH mass and a variety of host spheroid properties. Red points with error bars are the observations; red line shows the 
least-squares best- fit to each relation. Note that the quality of th e correlation between Mbh and concentration or Ser sic index observed is si gnificantly reduced 
relative to that in Graham et al. (2001); Graham & Driver (2006) when we adopt the Sersic index measurements from Kormendy et al. 12007). Black points and 
lines are the results of hydrodynamic simulations (see §|5J. The slopes of the simulated and observed relations are statistically identical in every case, with all 
normalization offsets owing to the small Faber-Jackson (M» (cr)) offset in FigurefT] 



outside the scope of this paper, as it would require a fully cos- 
mological prediction for merger rates (as a function of type, 
mass, and redshift) with high-resolution cosmological simu- 
lations such that the central galaxy structure and BH feedback 
effects could be modeled. The Mbh-M* and BHFP relations 
on which we focus, on the other hand, should be relatively 
robust to these effects (see also §|9). 

4.2. A Black Hole Fundamental Plane 

We wish to determine whether or not a simple one-to-one 
correlation between e.g. Mbh and a is a sufficient descrip- 
tion of the simulations, or whether there is evidence for addi- 
tional dependence on a second parameter such as R e or M* . 
The most efficient way to determine such a dependence is by 
looking for correlations between the residuals of the various 
projections of such a potential BHFP relation. Following Pa- 
per I, Figure [3] plots the correlation between BH mass Mbh 
and host bulge effective radius R e at fixed a. Specifically, we 



determine the residual with respect to the Mbh _ c relation by 
fitting Mbh(c) to an arbitrary log-polynomial 



(log(M BH )) = E [a„ logiaf] 



(2) 



allowing as many terms as the data favor (i.e. until A% 2 with 
respect to the fitted relation is < 1), and then taking 



Alog(M BH I cr) = log(M BH )- (log(M BH ))(er). 



(3) 



We determine the residual A \og(R e \ a) (or A log(R e \ a), for 
the stellar mass) in identical fashion, and plot the correlation 
between the two. We allow arbitrarily high terms in log(er) to 
avoid introducing bias by assuming e.g. a simple power-law 
correlation between Mbh and cr, but find in practice that such 
terms are not needed - as discussed in Paper I, there is no 
significant evidence for a log-quadratic (or higher-order) de- 
pendence of Mbh on a, R e , or M„, so allowing for these terms 
changes the residual best-fit solutions by \a. Of course, 
even this approach could in principle introduce a bias via our 
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TABLE 1 
BH-Host Correlations 



1 1 

Variables * 




Observed 






Simulated 




Nornializatiori - 


Slope^ 




Nornitiliztitioii 


Slope^ 


Sc titter 




o. jj zc w.v/u 


a (Y) -L n 10 fl 4/* 4- fl 1 Q 


21 


q i f. _|_ n ns 

O.IU w . Lf J 


? Qfl 4- IS 4- 1 1 

Z..y\i IH U.JO, U.JH ZT_ \J. 1 1 


21 


My, <r H 


O.Z4 ± U.Uu 


n </i j_ni*7 o 1 o i a co 
U.j4±U.l/, Z.ioitU.Jo 


U.ZZ 


/.yj ± u.uu 


U. 12. ± U. iz, 1 .4U ± U.4y 


n 1 q 

u. iy 




8.06 ±0.07 


1.78±0.40, -1.05±0.37 


0.25 


7.64 ±0.04 


1.50±0.22, -0.56±0.26 


0.21 


M* cr 2 


8.23 ±0.06 


0.71 ±0.06 


0.25 


7.92 ±0.04 


0.71 ±0.03 


0.21 


<T 


8.28 ±0.08 


3.96 ±0.39 


0.31 


8.04 ±0.06 


3.91 ±0.20 


0.31 


M* 


8.21 ±0.07 


0.98±0.10 


0.33 


7.85 ±0.05 


1.16 ±0.06 


0.23 


M dyn 


8.22±0.10 


1.05±0.13 


0.43 


8.18±0.06 


1.05 ±0.06 


0.28 


Re 


8.44±0.10 


1.33±0.25 


0.45 


8.48±0.12 


1.92±0.14 


0.56 


0c 








8.31 ±0.05 


1.77 ±0.08 


0.23 



For the variables (.v, y), a correlation of the form log(MBH) = a\og(x) + (3\og(y) + 8 is assumed, where the normalization is S and a, [3 are the logarithmic 
slopes.^ The normalization gives log(MBH/MQ) for a = 200kms~ 1 , M* = 10 11 Mq , Mayn = 10 11 Mq, R e = 5kpc, which roughly minimizes the covariance 
between fit parameters.' Errms quoted here for the BHFP relations in (cr, R e ), (M* , cr), and (M* ,R e ) include the covariance between the two slopes. 
Holding one of the two fixed and varying the other yields substantially smaller errors (typically ~ 5%). All quoted errors account for measurement 
errors in both Mbh and the relevant independent variables.^ The internal scatter is estimated from both the simulations and observations as that which 
yields a reduced x~ jv = 1 with respect to the given best-fit relation. Central potential c/> c , normalization at <p c = 10 6 km 2 s~ 2 . There are no observational 
measurements presently available to compare with this correlation. 

assumption of some functional form, and so we have also con- 
sidered a non-parametric approach where we take the mean 
(log(MBH)) in bins of log(cr). Our large set of simulations al- 
lows us to do this with very narrow binning, and we recover a 
nearly identical answer. 

The simulations show a highly significant correlation be- 
tween Mbh and R e at fixed cr, similar to the observed trend 
in residuals. We therefore introduce a FP-like relation of the 
form 

M BH <XCr Q /vf, (4) 

which can account for these dependencies. Formally, we de- 
termine the combination of (a, (3) which simultaneously min- 
imizes the x 2 / l/ °f the fit and the significance of the cor- 
relations between the residuals in a and Mbh (or R e and 
M bh)- This yields simila r results to the direct fitting method 
of iBernardi et alj d2003b) from the spheroid FP, which mini- 
mizes 

A 2 = [log(M Bll )-alog(a)-l3log(R e )-6]\ (5) 

and corresponds exactly to the method used in fitting the ob- 
servations in Paper I, to which we compare. This yields a 
best-fit BHFP relation 

log(M BH ) = 8.16 + 2.90(±0.38)log(cr/200kms- 1 ) (6) 
+0.54(±0.11)log(7v c /5kpc) 

from the simulations, similar to the observed BHFP relation, 
M BH oc a xo±lu fl0.43±o.i9 ( Paper Unsurprisingly, the slopes 
in the BHFP relation are close to those formally determined 
for the residuals in Figure [3] Figure [4] plots the residuals 
of Mbh with respect to these fundamental plane relations, 
at fixed R e and fixed cr. The introduction of a BHFP elimi- 
nates the strong systematic correlations between the residuals, 
yielding fiat errors as a function of a and R e . 

Given the definition of Md yn oc <r 2 R e , it is trivial to convert 
the best-fit BHFP relation in terms of a and R e to one in Md yn , 



obtaining 

log(M BH ) = 8.11+ 0.54 logCMdyn/TO 1 1 M ) (7) 
+1.82 logCcr/^OOkms- 1 ) 
= 8.03+ 1.45 log(M dyn /lO n M ) 
-0.91 log(/v,/5kpc). 

As is the case for the observed systems, the correlation with 
either a or R e at fixed Md yn is non-zero and highly significant 
(i.e. the BHFP is not simply a reflection of a simpler Mbh~ 
Md yn relation). 

At low redshift, cr, R e , and Md yn can be determined reliably, 
but at high redshift it is typically the stellar mass M* or lumi- 
nosity which is used to estimate Mbh- Therefore, it is interest- 
ing to examine the BHFP projections in terms of e.g. M* and 
cr or M* and R e . Repeating our analysis, we find in Figures [5] 
& [6] that the observations demand a FP relation over a simple 
MbhvM*) relation at high significance. We find a nearly iden- 
tical result using the dynamical mass, Md yn oca 2 R e , instead of 
stellar mass M*, as expected from the BHFP in terms of a and 
R e . The exact values of the best-fit coefficients of this BHFP 
determined from the observations are given (along with those 
of various other BHFP projections) in TableQ] 

Because the simulation merger remnant spheroids lie on a 
stellar-mass fundam ental plane very sim ilar to observed el- 
liptical galaxies (Robert son et alj l2006bl) . which tightly re- 
lates M*, R e , and cr, the BHFP in terms of any two of 
those variables can be easily converted to any other two. In 

other words, the two forms of the BHFP (Mbh oc a R l J 2 and 

9 1/2 

Mbh oc cr MJ ) are completely equivalent (the choice be- 
tween them is purely a matter of convenience), and it is re- 
dundant to search for a four-variable correlation (of the form 
Mbh oc a a R^ M*). As in the case of the observations (Paper 
I), we stress that no transformation (given the early-type FP 
relating these three variables) eliminates the dependence on 
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FIG. 3. — Upper: Correlation between the residuals in the Mbh _ <& relation 
and R e -cr relation, from our simulations (black points) and observed sample 
(red points with errors). At fixed a, systems with larger effective radii R e 
also have larger black hole masses Mbh- The fit to this residual correlation 
is shown with the black lines (± 1 a range in the best-fit correlation shown as 
dashed lines - note that they are strongly inconsistent with zero correlation), 
with the slope shown. The probability of the null hypothesis of no correla- 
tion in the residuals (i.e. no systematic dependence of Mbh on R e at fixed 
a) for the observed systems is shown (red P llu n) - the observations imply a 
secondary "fundamental plane"-type correlation at 3<x. Lower: Same, but 
considering the correlation between Mbh and a at fixed effective radius R e . 

two variables - i.e. no transformation from any one BHFP al- 
lows us to write Mbh as a pure function of either a, R e , M», 
orM dyn . 

5. IMPLICATIONS OF THE BHFP FOR LOCAL BH MASSES AND 
DEMOGRAPHICS 

Given that the "true" correlation between Mbh and host 
properties appears to follow a FP-like relation, it is natural to 
ask how adopting such a relation affects the estimation of BH 
masses from observed host properties. Figure|7]shows the ob- 
served and simulated systems in the fundamental plane. The 
relation appears to be a good predictor of Mbh over a large 
dynamic range, and there is no evidence for any curvature or 
higher-order terms in the relation (fitting e.g. a log-quadratic 
relation in this space yields A\ 2 < 1). As detailed in Table[T] 
the intrinsic scatter in the BHFP is small, typically ~ 0.2 dex, 
and in all cases smaller than the scatter in e.g. the M B h-<7 or 
Mbh~M* relat ions. 

However, as iNovak et al.l (|2006) note, minimizing the in- 
trinsic scatter does not necessarily maximize the observational 
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FIG. 4. — Correlation between the residuals in our BH "fundamental plane" 
relation Mbh °c cf 2 90 R® 54 and effective radius R e at fixed a (upper) or a at 
fixed R e (lower). Accounting for the joint dependence of Mbh on <r and R e 
removes the strong systematic dependencies in the residuals from Figure fJ] 
Cndl is large, meaning there is no further residual dependence). 



ability to predict BH masses. The BHFP relations depend on 
measuring two of either a, R e , or M*, and therefore introduce 
additional errors from the measurements of two (as opposed 
to just one) of these quantities. At low redshifts, it may be 
possible to obtain accurate measurements of both a and R e 
and therefore still obtain more accurate mass estimates (al- 
though we caution that several of the literature sources from 
which we compile observations differ by > 2a in some R e 
measurements, owing to various systematic issues such as the 
choice of observed bands). However, at high redshifts M* 
remains the most easily applicable proxy for Mbh, and it is 
not clear that the additional accuracy gained by introducing 
the R e term substantially improves the predictive power of the 
relation. 

Ultimately, the Mbh _ cr and Mbh _ M* relations are not 
much worse in a mean sense around M BH ~ 10 8 - 1O 9 M , 
with relatively small intrinsic scatter (see Table Q]). The rea- 
son these relations work as well as they do is that they are 
both nearly edge-on projections of the BHFP. Given a rela- 
tion Mbh oc ct 3 R° e 5 , it is not surprising that a is an acceptable 
proxy for Mbh in many situations (whereas the Mbh -Re cor- 
relation has quite large scatter, as R e enters with a relatively 
weak dependence in the BHFP). 

However, given that there is a systematic dependence on 
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FIG. 5. — As Figure[3] but considering the correlation between the residual 
in BH mass Mbh and stellar mass M* at fixed velocity dispersion cr (upper), 
or between A/bh and a at fixed M» (lower). Observed and simulated systems 
with larger velocity dispersions (deeper potential wells) at fixed stellar mass 
have more massive BHs, as do more massive systems at fixed a. A nearly 
identical result is obtained using dynamical mass M^ n instead of stellar mass 

Mi, : a residual dependence oc M° v f ±0M (P nu[[ = 0.0037) at fixed a and ot 



dyn 

= 0.0032) at fixed M iya , 



e.g. M„ or R e at fixed a which is captured only by the 
BHFP relations, we expect that the importance of estimat- 
ing BH masses from the BHFP will be enhanced at the ex- 
tremes of observed distributions. Figure [8] compares the lo- 
cations of outliers in the simulated and observed Mbh - cr and 
Mqh-M d yn relations with their locations on the BHFP rela- 
tions. Indeed, several systems which appear as outliers in one 
of the projections of the BHFP are no longer significant out- 
liers in the BHFP relation. This is true for a number of sys- 
tems on both the M BH -M* relation and the M BH - cr relation 5 . 
These systems typically have abnormally high or low velocity 
dispersions given their stellar mass, and therefore appear de- 
viant in the BHFP projections, just as such systems typically 
appear to deviate from the spheroid FP in projections such as 
the R e -M* or Faber-Jackson (M*-er) relations. Therefore, 
while the typical scatter about the mean relation is not dra- 
matically different for the BHFP (~ 0.2dex) compared to the 

5 NGC 1023, 3384, 4697, 5252, and Cygnus A have BH mass 
measurement errors of < 0.1 dex and measured masses which are 
0.32, 0.96, 0.81, 0.53 , and 0.34dex discrepant with the expectation from the 
ITremaine etal] f2002l) relation, but only 0.10, 0.70, 0.51, 0.44, and 0.23 dex 
discrepant with the BHFP expectation, respectively. 
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FIG. 6. — As Figuref4] but considering the residuals with respect the the 
BHFP defined in terms of stellar mass and <r. Placing the simulations and 
observations in the context of a FP relation eliminates the strong trends in 
their residuals. 



MftH-a relation (~ 0.3 dex), the tails of this distribution are 
substantially suppressed when we adopt the BHFP as a BH 
mass estimator. 

Given the potential importance of the BHFP for predict- 
ing the masses of BHs, especially in extreme systems, we 
should examine the implications that the relation has for 
the local demographics of BHs. There has been substan- 
tial debate recently about whether or not high-M* systems 
begin to deviate from the low-M* Faber-Jackson relation 
(see, e.g. iBovlan-Kolchin et all l2006t iBernardi et~aT] l2006t 
iBatcheldor et alj [20061) . If so, this implies that either the 
Mbh~M* or M Bll -cr relation must change slope at the high- 
est masses, with one of the two or some other relation being 
the conserved relati on. Because the d istribution of spheroid 
velocity dispersions ( Sheth et al. 2003) declines more steeply 
at high a than the galaxy mass or luminosity functions do at 
high M* (when BCGs are included), the assumption that the 
MflH-Mf, relation remains unchanged (the simplest expecta- 
tion from gas-free or "dry" mergers) naively predicts a much 
higher abundance of ver y high-mass (> 10 9 M(Ti) BHs (for an 
extended discussion, see lLauer et al . 2006b). 

Figure|9]compares the expected BH mass function (BHMF) 
from the observed distribution of spheroid velocity disper- 
sions and the Mbh _ cr relation with that expected from the 
early-type galaxy stellar mass function and the Mbu-M* rela- 




FIG. 7. — Masses of BHs in our simulations and from local measurements, compared to the expectation from the best-fit BHFP relations in a, R e and M,, a. 
The two agree well at all masses, without any evidence for curvature in the relations. The intrinsic scatter in Mbh at fixed a, R e or M» . a is estimated from the 
simulations to be ss 0.20dex (see Tableff}, which is consistent with the scatter in the observed points (given their measurement errors). 

the lCole et all (120011) or lJones etail (12006b /T-band luminos- 
ity functions from the 2dF and 6dF respect ively (converted 
to mass functions following Bell et al. (2003)), which include 
BCGs and extend to the low space densities of interest. (In 
any case we are simply attempting to highlight the key qual- 
itative behavior, that with the inclusion of BCGs the stellar 
mass function falls off more slowly than the velocity disper- 
sion function.) 

For now, we assume a one-to-one correlation between Mbh 
and a or M*. Given this, there would also be a one-to-one 
correlation between M» and a (i.e. er(M*) is given by match- 
ing the two distributions at fixed number density), so we can 
use the BHFP relation in the form Mbh oc M" to estimate 
the BHMF from this relation. Unsurprisingly, we find that 
the BHMF from Mbh _ cr cuts off rapidly compared to that 
from Mbh _ M» . With a mixed dependence on both a and 
M„, the BHFP relation predicts a somewhat intermediate case 
at high masses, although it is closer to the expectation from 
Mbh-M*. 

However, this treatment ignores the important fact that there 
is scatter in these correlations. To predict the BHMF from 
the distribution of velocity dispersions, we should properly 
convolve over the mean relation broadened by some (approx- 
imately lognormal) dispersion with width ~ 0.3 dex. The in- 
trinsic scatter is difficult to determine from the observations, 
if errors are not completely understood, so we adopt the in- 
trinsic scatter in each correlation estimated from our simu- 
lations (Table [TJ. Given this scatter, we re-calculate the ex- 
pected BHMFs, shown in Figure [9] Interestingly, the three 
predicted BHMFs are now almost identical, even at very high 
BH masses (> 5 x 10 9 ). This is expected - if we completely 
understand the correlation (and its scatter) between Mbh and 
either a or M* over the entire mass range of interest, then 
any projected version of the same BHFP must yield a simi- 
lar BHMF. The fact that the distribution of a cuts off more 
steeply than M» is compensated by the fact that the intrin- 
sic s catter in Mbh ~ o is s lightly larger than in Mbh-M* (see 
also [Marconi et al]|2004l) . This is roughly equivalent to the 
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FIG. 8. — Upper: Outliers in the local Mbh — cr relation from the observa- 
tions (red; defined as points > 3 a discrepant from the mean trend) and sim- 
ulations (black; defined as points 0.5 dex dis crepant from the mean trend) are 
plotted relative to the mean relation from Tremaine et al. (2002) (line) (left), 
and plotted relative to the expected BH mass from the fundamental plane rela- 
tion (right). Lower: Same, but for outliers from the local Mbh — Mi yn relation 
from Haring & Rix ( 2004). Most outliers in Mbh _ & or Mbh - 



-M. 



dyn are ex- 
plained by the BHFP relations we derive - i.e. they have abnormal values of 
cr or R e for their masses, but are not outliers in the BHFP relation. There are 
a small number of observed and simulated systems which are genuine out- 
liers from all relations - i.e. do have anomalous BH masses, but we note that 
there are no simulated or observed systems that are outliers from the BHFP 
relations and not also outliers from Mbh — a or Mbh ~^dyn- 



tion. We ad o pt the distribution of velocity dispersions a from 



ISheth et 



adopt tne c 
HJ (12001. 



and the e arly-type galaxy stellar mass 
function from iBell etail d2003l) with t he addition of BCGs 
at high masses from lLin & Mohrl (|2004) (yielding a shallower 
power law-like falloff at high masses, as opposed to an expo- 
nential cutoff). Note that a simi lar result is obtained estimat- 
ing the BCG mass function from lLauer et al.l (l2006bl) or using 
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FIG. 9. — BH mass function (cumulative number density above a given 
mass) expected from the obser ved early-type velocit y dispersion function 
(with the Mbh — o" relation from Tremaine et al. 1 2002)), the observed early- 
type stellar mass function (with the Mbh — M, relation herein, but this is quite 
similar to the Mbh ~Md yn relation from Haring & Rix 1 2004)), and from the 
joint distribution of M « and a (given the observed M» - a relations from 
IBernardi et al. 1 2003a) or implied by the two distribution functions, and our 
BHFP relation) or a and R e (with the observed Mbh _ Mj yll relation). The 
implied BHMF at high masses is very different for some of the various corre- 
lations, if we ignore the scatter in the relations (upper), with the BHFP repre- 
senting an intermediate case (albeit closer to Mbh — M* than Mbh _ c). How- 
ever, accounting for the intrinsic scatter in each correlation estimated from 
our simulations (lower) yields nearly identical BHMFs even at ~ 1()'°Mq. 
Larger scatter in the Mbh — <& relation at high masses reflects (and compen- 
sates for) the discrepancy at high-MflH ■ Dotted lines show the BHMF inferred 
from Mbh — c with a change of just 25% in the estimated intrinsic scatter 
(from a scatter of a = 0.31 to a = 0.27 , 0.36), which demonstrates that a 
small uncertainty in the intrinsic scatter (difficult to determine observation- 
ally) can dominate the choice of BH-host correlation adopted, even at very 
high M B h- 

statement that at large masses, where the relation between a 
and M* may change, the Mbh _ c relation must change cor- 
respondingly (following t he true BHFP re lation; for a more 
detailed comparison see lLauer et al. 2006b). Because the dis- 
tribution of a falls rapidly at high a, there is little contribution 
at Iow-Mbh from high-er systems, so a change in slope or in- 
crease in scatter in Mbh - cr both have the primary effect of 
increasing the expected number of high-mass BHs, reconcil- 
ing the BHMFs. 

The scatter is of critical importance at these masses: we 
consider the BHMF derived from Mbh — ct if we change the 
estimated intrinsic scatter by just 25% (i.e. within the range 
0.27- 0.36 dex), all within the range allowed by the present 



observations ( Trem aine et al.l [2002). and find that this rela- 
tively small difference in the intrinsic scatter estimate makes 
a larger difference at high Mbh than the choice of correlation 
(M B h - c, Mbh-M* or BHFP) adopted. This reinforces th e 
point empha s ized b y lYu"& Tremaine! (120021) ; lYu & Lul (120041) ; 
Tundo ^tail ((2006) that the estimated intrinsic scatter can 
dominate the demographics of high-mass BHs - accounting 
for this, the BHFP does not substantially change these esti- 
mates. 

6. THE PHYSICAL ORIGIN OF THE FUNDAMENTAL PLANE 

If BH growth terminates because of self -regulation, the fun- 
damental requirement is that sufficient energy be released to 
unbind the surrounding galactic gas. Given a radiative effi- 
ciency e r and feedback coupling efficiency rj, the energy cou- 
pled to the IGM from accretion onto the BH over its lifetime 
is simply 



E Bn = r)t r M BH c 



(8) 



and the binding energy of the gas in the center of the galaxy 
is 

E gas = 4>f^M,a 2 , (9) 

where is a constant that depends on the shape of the bulge 
profile (0=10.1 for a lHernquisi (1 19901) profile). In detail, the 
energy in Equation [S] is proportional to the accreted (as op- 
posed to the total) BH mass, and the binding energy of the 
bulge changes as a function of time during the merger. How- 
ever, this simple estimate is actually a reasonable approxima- 
tion to what occurs in the simulations, for two reasons. 

First, as is also demanded by numerous em pirical con- 
straints (e.g.. lSoltanlll982tlHopkins et alJl2007fl) . the major- 
ity (> 70-80%) of the BH growth occurs in the observed 
"qua sar" phase, w hich is constrained to be both short-lived 
(e.g. [M artini 2004) and near-Eddington (e.g. iKollmeier et al] 
2006), as occurs in the final growth phase near the end of 
the merger. Therefore, most of the final Mbh is accreted 
in the final e-folding of BH growth, over a Salpeter time 
fs ~ 4.2 x 10 7 yr. This is small compared to the cooling time 
of the galactic gas at ~ R e , so the approximation that the en- 
ergy released ~ Mbh is reasonable. Second, by the time of 
the quasar phase, the merging galaxies have coalesced, and 
the bulge is largely formed and in place, so the BH growth 
occurs in the relatively fixed potential of the remnant. 

Equating the energy needed to unbind the surrounding gas 
and terminate accretion yields the expected scaling 



M BH ~ 10 8 M G 



0.005 ( M * \( ° 
r)e r /f' UO n M / V200kms 



(10) 



where we adopt a iHernquisH ( 1 19901) profile and a typical 
e r = 0.1, rj ~ 0.05 (similar to the values adopted in our sim- 
ulations), which roughly reproduces the observed Mbh— M* 
normalization. Note that, in detail, the BH feedback need not 
be radiative - kinetic wind or jet feedback could inject com- 
parable energy (e.g ., Maraschi 20 041 ; rTavecchio et al.l 120041 : 
iMerloni et alj|2005l) . with some efficiency E = f] Mbh c 2 . This 
yields an identical Equation (TTOl . as 77 is functionally equiv- 
alent to the previous rje r - it does not matter in this deriva- 
tion whether the net feedback efficiency (77 e r ) represents a 
radiative or kinetic mode (or some sum of the two), or ex- 
actly what fraction of the final BH mass is accreted in a given 
growth phase, as the coupling efficiency simply serves to set 
the normalization of the MBH-host correlations. It should also 
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form 



M BH cx(M»a 2 ) c 



(11) 
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FIG. 10. — As Figuref5] but comparing the residuals in BH mass and host 
properties M*, a, and R e at fixed spheroid binding energy M*<t 2 . There 
is no significant evidence in the simulations or observations for a residual 
correlation in this space. This is because the fundamental plane in BH mass 
can be (approximately) represented as a "tilted" correlation between BH mass 
and bulge binding energy, Mbh oc (M* it 2 ) ' 71 . 



be noted that / gas changes by a large amount during the evo- 
lution of typical galaxies or (especially) mergers - the value 
relevant for Equation ( TTOb is some effective value in the cen- 
tral regions near the BH during the final e-folding(s) of growth 
(hence the notation fi as ), which does not necessarily trace the 
global or pre-active pliase gas fraction of the system (see §|7J- 
We therefore naively expect that the BH mass should scale 
with M»cr 2 . In Figure [10] we examine the residuals of the 
best-fit correlation between Mbh and this binding energy 
M* a 2 , in the manner of Figure [3] In this space, there does 
not appear to be any strong > 2 a evidence for a correlation 
of the residuals in Mbh(M* cr 2 ) with those of M*, cr, or R e . It 
seems that the correlation between BH mass and bulge bind- 
ing energy is in some sense more basic then the correlation 
between BH mass and e.g. M* or a. However, when we fit 
Mbh to a function of M* cr 2 , we do not recover Equation (TTOt 
- in fact, a linear proportionality between Mbh and M* cr 2 is 
ruled out at ~ 5 a in the observations (and > 10 ct in our sim- 
ulations). Instead, BH mass follows a "tilted" relation of the 



with a w 0.71. Furthermore, this relation is not exactly the 
same as the BHFP we recover from the observations, which is 
closer to M B h oc M°' 5 cr 2 . 

If we revisit our argument, we note that we have naively 
assumed that the accretion energy from the BH is coupled 
in an (effectively) infinitely short period of time and un- 
binds the surrounding gas. More properly, what occurs in 
our simulations is a pressure-driven outflow from the cen- 
tral re gions (e.g., Hopkins et al. 2006b; Hop kins & Her nquist 
2006). This implies that the necessary condition for self- 
regulation is the injection of sufficie nt momentum (at a rate 
p oc L/c; see e.g. iMurrav et al.ll2005l) to drive a galactic out- 
flow (total p oc M* cr) within the dynamical time near the ra- 
dius of influence of the BH (Rbh = GMbh/c 2 ), 



fdynC^BH) : 



^BH 

a 



'BH/ 

GM B 



(12) 



When the rate is below this threshold, it can drive material 
from the central regions where it is initially bound or infalling 
onto the BH, but the momentum coupled is insufficient to 
entrain the larger-scale material and the outflow fails to halt 
accretion (typical of the early-stage weak winds s een in our 
simula tions in earlier merger stages; see Figure 1 in lCox et"afl 
2006a). This condition gives us the requirement 



or 



pAt(x Ml H /a 3 oc M* cr, 



Mbh oc M l J 2 a 2 , 



(13) 



(14) 



similar to the observed BHFP. 

Ultimately, the qualitative conclusio ns of these d erivations 
are similar. Empirical constraints (e.g., Soltan 1982) demand 
that most BH mass is accumulated in high Eddington ratio 
quasar phases - so the details of accretion at lower Eddington 
ratios (whether set by e.g. the Bondi-Hoyle rate or other accre- 
tion mechanisms) are not important. This also implies that the 
BH mass is accumulated in a short period of time ~ 10 7 — 10 8 
years, such that whether growth is driven at the end of galaxy 
mergers, or in secular disk or bar instabilities, the environment 
(background potential) local to the BH is relatively fixed. The 
BH then self-regulates when it is sufficiently massive that its 
feedback energy can unbind infalling gas and halt accretion. 
Effects which deepen the local potential at the galactic center 
and increase the binding energy of gas near the BH will pre- 
vent gas from being unbound until the BH has grown more 
massive than it would otherwise. For example, as we show in 
detail in § [7] if spheroid progenitors are more gas-rich, there 
is more gas which can cool and form stars on small scales 
in the center of the spheroid, deepening the potential well 
there (i.e. yielding a more concentrated remnant with smaller 
R e and larger cr, while having little effect on the total stellar 
mass). This requires that any BH grow larger in order to un- 
bind nearby gas and halt accretion, in a manner similar to the 
scaling of Equation ( fl4b . 

We might also ask how sensitively this BHFP scaling or 
"tilt" is related to dissipational processes, as Ro bertson et al.1 
(2006b) demonstrate that the origin of the spheroid FP tilt 
lies essentially in the scale-dependence of dissipational pro- 
cesses such as gas cooling and star formation. From the 
deriva tion above, we would ex pect that although the processes 
from lRobertson etail d200 6b) might affect the structure of the 
merger remnants themselves, they should not change how, 
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fundamentally, the central BH self-regulates. Unfortunately, 
BH accretion is, itself, naturally a dissipational process, and 
therefore we cannot simply test this theory in our case by run- 
ning simula tions where gas dissipa tion is turned off. 

However, Rob ertson et al.l (|2006b) further show that it ac- 
tually requires substantial initial gas fractions for these dis- 
sipational effects to act - in mergers with very low initial 
gas fractions (< 20%) no "tilt" (in the spheroid FP) is in- 
duced, consistent with requireme nts from the observed phase 
space densities of ellipticals (e.g.. lHernquist et aT1ll993h . We 
therefore briefly consider just a set of simulations with initial 
/gas = 0.05. This is sufficiently low that the remnants act dis- 
sipationlessly, and lie on the v irial relation as opposed to the 
spheroid FP (see Figure 10 of iRobertson et al.l d2006b)). but 
sufficiently large that we do not need to worry about artifi- 
cially "strangling" the BH by giving it insufficient gas to ac- 
crete (given typical Mbh ~ 0.001 M*, the BH need only have 
access to ~ 2% of this gas to grow normally). 

We find that these merger remnants obey a similar BHFP re- 
lation to their high-/ gas counterparts, implying that so long as 
feedback-driven self -regulation (as opposed to e.g. gas starva- 
tion) determines the final BH mass, these scalings are robust. 
However, the significance of the preference for e.g. a BHFP 
relation as opposed to a simpler Mbh oc M* or Mbh oc a A re- 
lation is greatly reduced. This is because, without significant 
effects of dissipation to change the central phase-space struc- 
ture and potential depth of the remnant, the velocity disper- 
sion a and effective radius R e are simply set by the violent 
relaxation of the scattered stellar disks. The general scalings 
of the BHFP are not, then, unique to gas-rich progenitors, but 
their significance and the importance of accounting for the 
observed dependencies are so. 

Indeed, regardless of prescriptions for BH accretion, feed- 
back, and star formation, the above derivations (and qual- 
itative dependence on galaxy properties) depend on only 
three relatively robust assumptions - that BH growth is dom- 
inated by bright, high-Eddington ratio phases, that feed- 
back from accretion affects the gas on small scales around 
the BH (sensitive to the central potential of the galaxy), 
and that some sort of heating, momentum coupling, or 
effective pressure eventually halts accretion. This does 
place some constraints on scenarios for BH growth, how- 
ever. In models where, for example, the BH accretion 
rate is a pure function o f the galactic star formation rate 
(see e.g.. iKawakatu etail l2006t iMonaco & Fontanotl l200l 
iGranato et all2004t ICattaneo et al.ll2.Q05l although in some of 
these cases the Eddington limit is still preserved as an up- 
per limit), it is difficult to explain how the BH mass would 
be sensitive to the central potential in the manner of the ob- 
served BHFP, and not simply trace the galaxy stellar mass (i.e. 
a yield pure Mbh — M» relation, which the observations disfa- 
vor at ~ 3fi). 

7. DRIVING SYSTEMS ALONG THE FUNDAMENTAL PLANE 

We have thus far considered systems in terms of the observ- 
able properties of the remnants, specifically quantities like 
M», er, and R e . We now turn to the "theorist's question" - 
namely, how are the positions of systems on the Mbh~M*, 
Mbh _ c, and BHFP relations affected by theoretical quanti- 
ties or initial conditions? 

Figure QT| considers the residuals in the Mbh _ M* and 
BHFP relation for merger remnants with different virial ve- 
locities, redshifts, and orbital parameters. Note that "redshift" 
in this context simply refers to the characteristics of the pro- 




Orbilal Parameters 



Orbital Parameters 



FIG. 1 1 . — Residuals of the Mbh -Af» (left) and BHFP (right) relations as 
a function of various initial conditions or theoretical quantities. V c is the halo 
virial velocity, z is the redshift to which the progenitor disks are initialized 
(i.e. their concentrations and scale lengths are rescaled by a small amount to 
match disks at these redshifts), and the orbital parameters shown span a small 
set ranked from those that yield the least (0) to most (3) rotation in the merger 
remnant s (average (V e /cr)* = 0.06, 0.25, 0.41, 1.03, respectively; for details, 
see Cox et al. 2006b). By changing e.g. the circular velocities, initial disk 
formation redshifts, or orbital parameters of simulated mergers, we can drive 
changes in the Mbh - A/* relation. However, these can all be understood in 
terms of how they change the galaxy structure - i.e. they drive changes in a 
at fixed M* , while preserving the BHFP relation. Higher angular momentum 
mergers also produce more disk-like remnants - but when we are careful to 
exclude the contribution from rotation to the velocity dispersion cr, the objects 
lie on the expected BHFP. 



genitor disks (which are initialized to resemble disks at low 
redshifts z ~ 0, or higher, z = 2, 3, 6). We select each set of 
simulations in which all parameters are identical except that 
plotted in the figure, then consider the residual with respect to 
the mean relation for just those simulations (each set of points 
of different color plotted in the figure represents one such set 
of simulations). 

In terms of Mbh/M*, there are weak (~ 0.2-0.3dex over 
the maximal range spanned by the simulations) trends towards 
lower Mbh/M* for higher redshift and larger angular momen- 
tum mergers, but the systems lie on the same BHFP regard- 
less of V v i r , z, orbital parameters, or any other initial quanti- 
ties we vary. The trends in Mbh/M. simply reflect changes in 
the structural properties of the remnants - for example, par- 
ticular orbital parameters produce smaller bulges and lower 
values of a at fixed M». This is expected - because most of 
the BH growth occurs over a short period of time at the end 
of a merger or other phase of activity (the last 1-2 e-folding 
times) during which the environment is relatively dynamically 
settled, changing initial conditions should only affect Mbh in- 
directly by altering the central structure of the remnant. The 
lack of a strong trend in y v ; r simply reflects the fact that the 
central regions which set the potential depth of relevance for 
determining Mbh are very much baryon-dominated. This is 
not to say that M B h does not scale with V v ; r , in a mean sense, 
but simply that the dependence on V v ; r is subsumed in the 
more direct dependence on M % and a, which themselves can 
depend on halo mass or V v j r . 



Modeling the Black Hole Fundamental Plane 



13 



1/2 9 

For this reason, at fixed MJ a , we expect BH mass to be 
independent of halo mass and redshift. Of course, M* and a 
will, in the mean, scale with Mh a i , which could yield evo- 
lution in the M BH - M^. d \ relation. For a bulge-dominated 
system (i.e. ignoring the complication that the same mass 
halo could host a disk-dominated system with a much smaller 
BH), M* = /*(fii,/Q m )A/haio, where /*(Mhaio,z) is the typi- 
cal fraction of baryons incorporated into the galaxy. In the 
simulations, we find a rough correlation a cx v max (similarly, 
v c oc v max for our progenitor disks, although that is by con- 
struction), where v max is the maximum halo circular veloc- 
ity (cx Vyh-c 1 / 2 , where c(M na i ,z) is the halo concentration), 
modulo the effects of e.g. gas fraction and orbital param- 
eters changing a at fixed Mh a i - Observationally, a similar 
mean correlation (a w 0.66 v max , nearly id entical to the best- 
fit normalization in our s imulations; e.g. iRronawitter et al.l 
120001: iGerhard et all 1200 ll) is found. Given c cx M^ a ° lo 13 (l + 
z)~ l dBullock etal.ll2001l) . this implies M BH cx a(z) f^ 5 M l h ^, 
where a(z) represents the weak remaining redshift evolution 
term, a = (Vt m A c (z)/fl m (z) 18-7T 2 ) 1 / 3 , which changes by only 
~ 20% from z = 0-6. We therefore expect that evolution in 
the Mbh-M^io relation will be dominated by the effects of 
evolution in typical gas fractions and remnant structural prop- 
erties on a (changing a at fixed Mh a i in a systematic sense) 
as well as cosmological evolution in typical baryon incorpora- 
tion fractions /* and/or bulge-to-disk ratios in galaxies hosted 
by halos of a given mass. To the extent that such effects occur, 
they of course should be also traced in some mean evolution 
in e.g. the a — v max relations. 

We have also studied a large subset of simulations of dif- 
ferent masses and gas fractions with a variety of prescriptions 
for feedback and winds from star formation. These results 
will be discussed in detail in Cox et al. 2007 (in preparation), 
but we briefly summarize their relevant conclusions. The in- 
clusion of massive stellar winds (with high mass loading ef- 
ficiencies rjy, > 1, where M w i n(J = -q^M*) can affect the struc- 
ture of merger remnants, although not necessarily in a mono- 
tonic or easily predictable fashion. For example, strong winds 
can remove gas from the central regions of the galaxy and 
yield a lower effective gas fraction / gas , but they can also cy- 
cle such gas at earlier stages, preventing it from immediately 
being turned into stars before the final merger, and actually 
raising the effective / gas at the final coalescence. In any case, 
the effects are usually small (especially in the most massive 
galaxies of interest here) or comparable to those owing to the 
choice of orbital parameters. 

More important, regardless of the stellar feedback effi- 
ciency, star formation alone cannot prevent gas from accreting 
onto the central BH (especially given that only ~ 0.1% of the 
galaxy mass in gas needs to reach the BH to affect the Mbh- 
host correlations), and in fact strong stellar winds can (either 
by cycling gas as above, or by shocking and increasing the 
supply of low angular momentum material) greatly increase 
the fuel supply for accretion at the galactic center. The cor- 
relation between BH mass and host properties is therefore, 
regardless of the prescription for stellar feedback, still set by 
the local self-regulation of the BH, and obeys (in all our sim- 
ulations) an identical BHFP relation. Stellar winds can, in 
principle, influence the final BH mass, but only indirectly by 
affecting the structure of the remnant, in the same manner as 
changing progenitor structural properties, orbital parameters, 
and gas fractions. 

Figure [12] considers the trend in Mbh / M» as a function 
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FIG. 12. — Trend in M^/M* as a function of simulation initial disk gas 
fraction. Upper panel plots Mbh/M* for three suites of simulations (solid, 
dashed, and dot-dashed lines). Each of the three is a set of otherwise iden- 
tical mergers of Milky Way-like systems, varying only the gas fraction with 
values /g as = 0.05, 0.2, 0.4, 0.6, 0.8, 1.0. The three suites consider three 
different orbital configurations for the merger. Lower panel considers t hese 
residuals and those from our other simulations versus / gas as in Figure [TT] 
Black line is a log-linear fit (dashed show ± 1 cr), green line the more accu- 
rate fit log(M BH /M*) = -3.27+0.36erf[(/ gas -0.4)/0.28]. High-/ gas mergers 
produce much larger BHs than low-/ gas systems. 

of the initial gas fraction of our simulations. In contrast 
to trends with V v h, orbital parameters, or the evolution of 
disk structural parameters with redshift, the dependence of 
Mbh/M* on / gas is quite strong, varying by nearly an or- 
der of magnitude from low (/ gas < 0.2) to high (/ gas > 0.8) 
initial gas fractions. We consider this in detail by exam- 
ining a small case study set of simulations. We construct 
a fiducial set of simulations of Milky-Way like initial disks 
(Kir = 160kms _1 ), and collide them in otherwise identical 
mergers except for varying the initial disk gas fractions with 
values / gas = 0.05, 0.2, 0.4, 0.6, 0.8, 1.0. We construct three 
such suites, each with a different orbit (roughly bracketing the 
extremes of possible merger configurations). Figure [121 also 
shows the trend of Mbh/M* in these simulations - it is clear 
that, all else being equal, larger values of / gas drive the sys- 
tems to larger Mbh/M*. 

However, the trend here does not resemble the simple 
Mbh oc / gas scaling that we naively predicted in Equation ( TTOb 
by demanding that the BH be able to unbind the entire initial 
gas content of the galaxy. In fact, Figure[T3]shows the correla- 
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FIG. 13. — Left: Correlation of Mbh with M* in simulations with different gas fractions / gas = 0. 2, 0. 4, 0.8, 1 .0 (orange, green, blue, and black, respectively; 
upper), and residuals of Mbh — Af» as a function of / gas for our case study simulations from Figure[T2](ZoM.'er). Center: Same, but for the A/bh — o correlation. 
Right: Same, but for the BHFP relation. Varyin g / ga s does not move the remnants off the BHFP relation, but does systematically shift them with respect to M 4 
and (j. The trend of Mbh/M* with / gas in Figure [l2l must therefore fundamentally relate to how / gas modifies structural properties like a (preserving the BHFP), 
not to the naive expectation that larger / gas translates to "more material to be unbound." 

tion between Mbh and M* a 1 for simulations of different gas 
fractions, and there is no systematic trend with / gas . Likewise, 
the remnants lie on the BHFP regardless of their gas fractions 
- i.e. the change in Mbh/M* can be entirely accounted for by 
the change in a at fixed M*. This should, perhaps, not be sur- 
prising - our earlier derivation neglected the fact that, by the 
quasar phase and epoch of final BH growth, the large major- 
ity (> 80-90%) of / gas has already been turned into stars or 
ejected by stellar winds. The trend in M BH /M* with / gas must 
therefore be primarily driven by how / gas changes the struc- 
tural properties of the remnant, a more subtle effect than the 
naive "amount of material to be unbound" expectation. 

In order to understand this strong dependence of Mbh/M* 
on /g as , we consider the structural properties of the merger 
remnants in our case study in Figure [14] On inspection, it is 
clear that the trend in Mbh/M» is one of increasing Mbh at 
fixed M*, as M* is nearly constant with / gas (there is a weak 
trend, as not all of / gas is converted to stars, but this changes 
M* by < 0. 1 dex from / gas = to / gas = 1 ). However, there is a 
strong trend in R e and a with / gas , as increasing the amount of 
dissipation (through / gas ) yields more concentrated remnants. 

We can understand this behavior with a simple toy model. 
Assume that the stars formed in the disk(s) before and dur- 
ing the merger (i.e. those that will a ct dissipationlessl y) are 
scattered into a typical bulge with a Hernquist ( 19901) pro- 
file with scale length /? c (/ gas = 0), independent of the cen- 
tral gas content, and that a fraction fj, of the initial gas mass 
(M gas = /gasM ga i) survives to the late stages of the merger and, 
via dissipation, falls to the center of the galactic potential. 
There, it will form a highly concentrated central stellar com- 
ponent (scale length <C ^ c (/ gas = 0); for simplicity we take it 
to be effectively a point concentration). The total enclosed 
mass as a function of radius is then 




FIG. 14. — The structural properties o f the merger remnants from our set 
of gas fraction case studies from Figure [T2] Changing / gas has almost no 
effect on M», but by increasing the amount of dissipation and fraction of 
stellar material formed in the final, central starburst, increasing / gas produces 
more concentrated remnants with smaller effective radii R e and larger central 
velocity dispersions a. Red lines show the expectation of a toy model in 
which a fraction ~ 0.5 (s olid ; ±0.05, dashed) of the gas participates in a 
central starburst (Equation ll6t . 

This yields a half-mass effective radius R e of 
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Figure [14] compares this simple expectation for R e and a 2 oc 
M*/R e with that from the simulations as a function of / gas 
- for a representative /i = 0.5, our toy model describes the 
simulations quite well (until / gas — > 1, where our assump- 
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tion that the inner stellar component is infinitely concentrated 
breaks down; this "saturation" in a reflects the fact that even 
extremely gas-rich systems will still turn much of that gas 
into stars that scatter into some large orbits, and may be im- 
portant for e.g. the steep observed cutoff in the observed a 
distribution). Note that this scaling is not very sensitive to 
our assumption about the exact profile shape - assuming the 
bulge follows an exact r 1 / 4 law or adopting a different in- 
ner power-law slope changes the predicted scaling by only 
~ 10-20%. Given this change in R e at fixed M* owing to in- 
creasing / gas , the fundamental plane implies that Mbh should 
increase, roughly as ~ R~ l . Figure Q3] plots the dependence 
of Mbh/M* on / gas , compared with the expectation from this 
simple model. From the agreement here, and the fact that Fig- 
ure Q3] finds no change in the BHFP with / gas , we conclude 
that more gas-rich mergers drive evolution in Mbh by produc- 
ing more concentrated remnants with smaller R e and larger a 
at fixed M„, and therefore larger M*£r 2 . 

There is one important caveat to our discussion of disk 
gas fractions / gas . Lacking a full cosmological simulation in 
which to determine how gas continuously accretes onto the 
disks, we have simply referred to / gas as the initial gas frac- 
tion in our simulations. Of course, during a simulation, / gas 
will decrease as gas is turned into stars, so that the actual gas 
fractions by the time the systems merge may be substantially 
lower than the numbers we quote. In Figure[l5] we reproduce 
our plot from Figure[T2]of the residuals in Mbh/M. as a func- 
tion of the initial gas fraction of the simulations. However, 
we also return to the simulations and measure the gas fraction 
for each at a set of uniform times Af = 1.0, 0.5, and 0.2 Gyr 
before the final merger (defined for convenience as the coa- 
lescence of the two BHs). As the simulations approach the 
final merger event, it is clear that the trend of residuals with 
/gas is qualitatively unchanged - however, the absolute values 
of /gas systematically decrease. By Af = 0.2 Gyr, the gas frac- 
tions are systematically lower by a factor ~ 2— 3. Therefore, 
the exact values of / gas which we quote should not be taken 
too literally - if gas is accreted in the real universe such that 
/gas changes less rapidly in earlier stages of a merger, then the 
initial gas fraction need only be as large as / gas ~ 0.3 to be 
equivalent to our most extreme / gas = 1 cases (with our more 
typical /gas = 0.4-0.5 cases corresponding to rather moderate 
pre-merger or Af = 0.2 Gyr gas fractions of ~ 0.2). 

8. IMPLICATIONS FOR THE REDSHIFT EVOLUTION IN BH-HOST 
RELATIONS 

8.1. Empirical Predictions 

Given that the BHFP appears to be robust against all var- 
ied quantities in our simulations, we expect that it should be 
preserved at all redshifts. However, this implies that, at fixed 
M*, evolution with redshift in the typical velocity dispersions 
and/or effective radii of spheroids will also manifest as evolu- 
tion in the typical Mbh/M* rel ation. 

It is observed empi rically dTruiillo et alj 120061) and ex- 
pected theoretically dKhochfar & Silkl 120061) that high- 
redshift spheroids will be more compact at a given stel- 
lar mass M* t han th eir low-redshift analogues. Specifically, 
Truiill oet al.l (H006) compile a number of measurements of 
the evolution, relative to z = 0, of the effective radii of 
spheroids (defined as systems with Sersic indices n s > 2.5) 
at fixed stellar mass M* > 3 x 10 10 M o and > 6 x 1O 1O M , 
corresponding to typical L* galaxies at most redshifts. If the 
BHFP is preserved, this necessarily implies evolution in the 
M BH -M* relation, of the form Mbh/M* oc R e {M*)~ l (see Ta- 



bleQ]for exact values). 

Figure [16] plots this e xpected evolution in Mbh/M. from 
the iTruiillo et alj ((2006) measurements, nor malized to the 
value observed at the same stellar mass by lHaring & Rixl 
(2004) atz«0. For comp arison, w e plot t he estimated evo- 
lution in Mbh/M* from Pen g et all (120061) (specifically, we 
adopt their early-type template to convert their measured lu- 
minosities to s t ellar m asses), and limits on this evolution from 
Hopkins et al. (2006d), normalized to the same local value. 

We can also attempt to empirically infer the evolution in 
Mbh/M* by considering the clustering of quasars as a func- 
ti on of redshift. Essentia l ly, thi s expands on the measurement 
in lAdelberger & Steidell (120051) . A more deta iled discussion 
of this ( and t he samples we co nsider) is given in Hopki nTet alj 
(2007e|) and lFine et al.l ([2006), but we briefly review it here. 

We consider a quasar sample (typically observed near ~ L* 
in the quasar luminosity function) for which BH masses or 
typical Eddington ratios have been directly measured (for 2dF 
quas ars, we directly a dopt the observed BH mass distributions 
from lFine et alJ (l2006h . otherwise we use tho se determined as 
a func tion of quasar luminosity and redshift in Kollmeier et al. 
(2006)). We then use the observed clustering properties of 
that sample to infer a characteristic host halo mass - in other 
words, match the observed large-scale bias of the quasar pop- 
ulation to the average (number-density weighted) large-scale 
bias for halos in some mass range (at the same redshift). We 
calculate the ex pected large s c ale bi as of halos as a function of 
mass following Mo & White! fl 996) with the improved fitting 
formulae from Sheth et al If we assume the ratio of 

galaxy to halo mass for ~ L* galaxies does not evolve much 
with redshift (which appe ars to be obseryationally confirmed 
to at l east z ~ 1; see e.g . iHevmans et alj|200"6l IConrov et alj 
120071: IZheng et~a l. 2007), then we have (implicitly) obtained 
the average host galaxy mass of these quasars, for which we 
know Mbh - i.e. an estimate of the mean Mbh/M*(z). Of 
course, the bias as a function of halo mass will depend on 
the cosmology adopted (specifically the value of cr%); this and 
e.g. the BH mass measurements from quasar spectral energy 
distributions all introduce fairly large systematic uncertainties 
(at least factor ~ 2) in the absolute implied value of Mbh /M* . 
However, to the extent that we are interested only in the rel- 
ative evolution of this with redshift, these uncertainties are 
much smaller. 

Compiling a number of measurements of quasar clustering 
as a function of redshift, the inferred evolution in Mbh/M. 
is shown in Figure [17] This provides a c ompletely indepen - 
dent measurement of Mbh/M* from that o flPeng et al.l (2006), 
with entirely different systematics, but nevertheless is in rea- 
sonable agreement with their estimates, and with our simple 
expectation from the BHFP relation. Of course, the present 
observations are not sufficiently robust to distinguish between 
evolution by a factor of ~ 2 versus evolution by a factor of 
~ 3 at high redshifts, but different probes seem to suggest a 
roughly comparable effect to that predicted. 

8.2. A Dissipation-Driven Explanation 

Given the BHFP, we have the empirical expectation that, at 
high redshift as spheroids become more concentrated, Mbh 
must be larger at fixed M* . However, this does not explain 
what physically drives these trends. In § [7] we showed that 
increasing the gas fractions of merger progenitors has both of 
these effects: namely, that by increasing the amount of dissi- 
pation, more centrally concentrated remnants with smaller R e , 
higher a and larger Mbh at fixed M* are produced. It is both 
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FIG. 15. — Dependence of M%n/M* on / ga s, as in Figure [T2l Panels plot this as a function of / gas measured at the initial time in each simulation (left) and at 
times At = 1.0, 0.5, and 0.2 Gyr before the final merger and q uasa r phase. Black lines in each plot the best-fit trend (solid; with ±lcr range in dashed lines), and 
red lines plot the expectation from our simple model (Equation [16), assuming that the BHFP is always preserved and that a fraction fi of the gas mass measured 
at each time will participate in the final, central starburst (/1 is smaller at early times because much of this gas will form stars in the two disks well before they 
merge). Regardless of when / gas is defined, the trend is similar, and evolution in Mbh/W* is driven by preserving the BHFP and increasing a at fixed M* in the 
same manner. 




FIG. 16. — Evolution in Mbh/W* as a function of redshift. Red stars plot the evolutio n expected at fixed M « (filled ~ 6 X 10 10 Mq; open ~ 3 X 10 10 MQ)given 
the observed evolution in the effective radii R e (M t ) for systems of this stellar mass from Trujillo et al. ( 2006) and our fundamental plane relations (Table[T). Other 
point s compare observational estimates of this ev olution: black s quares are direct measurements of high-redshift BH masses and host luminosities <Peng et al] 
2006), blue circles are upper limit estimates from Hopkins et al. (2()06d) based on observed spheroid and BH mass functions. All points are normalized to the 
same local value (black star) from Haring & Rix (2004). Solid and dashed lines show the predicted evolution from an a priori model in which the evolution 
in R e (Mt) (and, correspondingly, Mbh(^*)) is driven by increasing disk gas fractions with redshift, given the best-fit scalings of AfBH(/gas \M r ) (solid; dashed 
show ±lcr). The BHFP predicts that as high-redshift spheroids are more compact, Mbh/M* must rise, in good agreement with the observations; the trend is 
driven by more gas-rich progenitors in spheroid-producing mergers. To contrast, dotted line shows the expectation of the simplified (but common) semi-analytic 
assumption that Mbh oc a 4 oc V c 4 (assuming M„ (A/ V j r ) does not change much with z for ~ L„ galaxies). 



expected and observed that high-redshift disks are character- 
istically more gas-rich, as star formation has simply had less 
time to operate - this provides a potential a priori physical 
motivation for the evolution we saw in § 18.11 

We therefore briefly examine how systematic evolution in 
the gas fractions of the progenitors of gas-rich, spheroid form- 
ing mergers (i.e. disk galaxies) should affect the resulting 
Mbh -host galaxy correlations of the merger remnants (by, for 
example, increasing the amount of dissipation and therefore 



forming more compact, higher-cr remnants). To begin, we 
need to construct an estimate for how the gas fractions of 
typical disks evolve as a function of redshift. Fortunately, 
the structural properties (i.e. basic kinematics and the Tully- 
Fisher relation) of disks appear to evolve relatively weakly 
with redshift, so we can at least attempt to simply scale 
up from the properties of local disks (using e.g. their star- 
formation histories). Traditionally, the star formation histo- 
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FIG. 17. — As Figure [T6l (solid and dashed lines are our prediction and 
±1<t range, respectively; dotted line is an alternative Mbh °c V£ model), 
but the observationally estimated evolution of A/bh/M* (points) is derived 
from meas urements of quasar clustering as a function of redshift, following 
IFine et al] 420061) . For quasars of a given (measured) BH mass, their cluster- 
ing implies a characteristic host halo mass - if the ratio of that halo mass to 
stellar mass d oes not change with red shif t, then the implied M^u/M* evolves 
as plotted (see Hopkins et al. ( 2007e) and Fine et al. ( 2006) for details). Mea- 
sur ements of q uasar clustering are taken from Fine et al. 1 20061 red squares) 
and Porciani & Norberg (2006, orange stars) from the 2dF, and Myers et al. 
120061 blue circles) from the SDSS. 

ries of local disks are fitted to r-models, of the form M* <x 
exp[-(f-?,)/V], where r is some characteristic timescale and 
f, is an initial time of formation. This is, of course, a non- 
unique parameterization of the star formation history, but nev- 
ertheless appears to be a reasonable descrip tion of average 
stellar populations (e.g., Noe ske et al.l f2007). For a system 
with stellar mass M*(z = 0) = M*(f = f H ) Oh being the Hubble 
time at z = 0), this implies a normalization 

M*W = M * (Z = 0) (l-exptO-feVT])" 1 exp[(f,-f)/T]. 

(17) 

It is also well-establi shed that disks o bey a Kennicutt-Schmidt 
star formation law dKennicutC l 998 ) relating the surface den- 
sity of star formation to the gas surface density as Esf ~ ^g a t- 

The total SFR M* oc E SF /^, and £ gas cx M g . ds R~f, where 

-^gas = /gas^tot = /gas (1- /gas) -1 A^* and R d is the disk scale 
length. Knowing, therefore, how the mean star formation rate 
in disks of a given mass evolves, and how this scales with 
the surface density of gas, we need only some estimate of 
their characteristic sizes (Rd) to infer the (mean) evolution 
in the total gas mass content. If the baryonic Tully-Fisher 
relation do es not evolve with redshift, as suggeste d by ob- 
servations (tConselice et al . 2005; Flores et al. 2006) - i.e. the 
kinematic structure of disks does not strongly evolve, then 



we have M u 



l disk 



(Bell & de Jong 2001; McGaugh 2005, 



Vdi s k being the disk circular velocity = y/ GM lot /R d ) and there- 
fore M tot cx R d . Combining the Tully-Fisher relation with the 
Schmidt-Kennicutt law then implies 
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If we demand that the normalization of the Schmidt-Kennicutt 
law agree with the normalization from the fitted r model, we 
then arrive at an equation for the implied evolution in / gas , 
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/o 1 - 4 fl- exp[ft-f H )/T] 
-exp[(f,-f)/r] 



l-/o V 1 



exp[(f H -f)/r], (19) 
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FIG. 18. — Upper: Comparing the predicted evolution in characteristic 
disk gas fractions (as a function of stellar mass) from our simple toy model 
(in which we integrate backwards in time the best-fit r-model star forma- 
tion histories of local disks) with that observed. We compare the expected 
mean / gas (M t ) with estimates at z = (blue points; Bell & de Jondl200H 
Kannappan 2004; McGaugh 2005, as diamonds, squares, & , circles respec- 
tively ), z ~ 1 (green squares ; estimated from Ha luminosities Sha plev et all 



2005), and z ~ 2 (red circles Erb et al. 2006). Lower: Evolution in / gas from 
our simple empirical model for a constant M, = 1O"M0, compared to a 
simple exponential history with e-folding time f gas (as labeled). 



where /o = / gas (z = 0). Adopting the measured best-fit r and 
z = gas fraction fp as a f uncti on of z = stellar mass M* 
from iBell & de Jong (l2000h and iKannappanl d2004l) . respec- 
tively, we finally obtain an expected typical / gas and SFR as a 
function of disk stellar mass at any cosmic time t . 

Figure [18] compares the expected f gSLS (M^(t)) from this pa- 
rameterization with that observed at a number of different red- 
shifts. We also compare this estimate to another, even sim- 
pler parameterization we could adopt: assuming an isolated 
disk obeying a r-model, with / gas = 1 at t = and / gas sa 0.1 
(appropriate for a Milky Way-like ~ disk today) at z = 0. 
This implies an exponential growth of / gas with lookback 
time, with e-folding time ~ 6 Gyr. The results are similar, 
and appear to describe the observed / gas evolution reason- 
ably well. We have also checked that the expectation from 
Equation dT8l is consistent with observed specific star for- 
mation rate s as a function of M* from z = 0-3 (see, e.g . 
iBauer et alj|2005l iFeulner et al.ll200H iPapovich et aT][2006l) . 
and find good agreement even up to z ~ 3 (which should not 
be surprising, as this essentially just says that the r model 
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FIG. 19. — Upper: Faber-Jackson (M„ —a) and R e -M» relations in our 
simulations as a function of gas fraction (/ gas = 0.2, 0.4, 0.8, 1.0 in orange, 
green, blue, and black points, respectively, with the best-fit trend to each as 
the dotted line of corresponding color) . At fixed M» , higher-/ gas systems have 
smaller R e and larger a, as in Figure [14] Lower: Given our empirical esti- 
mate of the evolution in / gas at fixed M* = 10" Mq, and the dependence of 
R e (M t ) on /g as (upper panel & Figure [Ml . solid line plots the expected evo- 
lution in R e (M t ) with redshift. Points compare the observed evolution from 
ITruiilloetalJ f2006), for M, > 3 X 1O 1O M (filled) and M, > 6 X 10 I0 M Q 
(open). The increasing gas content in high-redshift merger progenitors pre- 
dicts their size evolution at fixed stellar mass. 



is indeed a reasonable description of the mean star formation 
history). 

Given this evolution in / gas , we are now in a position to es- 
timate how this should change the effective radii and velocity 
dispersions of merger remnants, and (given the BHFP) the av- 
erage value of Mbh at fixed stellar mass M*. Figure [T9lshows 
the evolution of the M» — a and Rg-M* relations with / gas , 
from our simulations, where the behavior is similar to that we 
discussed in §|7] Combining the trend in R e (M*) as a function 



of /g as from our simulations, and the trend in / gas (z) above, 
we also then predict how R e (z) s hould evolve, at fix ed M*. 
Comparing this to the results from Truii llo et al.l (|2006) shows 
excellent agreement given our simple estimate of / gas (z). 

Likewise, either directly adopting our fitted trend of 
Mbh/M* as a function of / gas from Figure [T2] or using our 
estimate of how R e (and, correspondingly, a) evolve at fixed 
M* with / gas and applying them to the BHFP relation Mbh oc 
M° 7 c 15 oc M^R' 1 , this predicts that the mean M B h/M* 
should increase with redshift. Figures Q~6] & [T7] show the 
evolution in Mbh/M» expected from this simple derivation. 
Again, the agreement with the observationally estimated rate 
of evolution in Mbh/M„ is good. 

We do caution that our estimate of the mean evolution in 
/gas, while consistent with the observational constraints on gas 
fractions and specific star formation rates at all redshifts with 
which we compare (z < 3), is only intended as a rough lowest- 
order approximation to realistic disk evolution. It is of course 
possible that other disk properties, such as the Tully-Fisher 
relatio n, begin to evolve at high redshift (see e.g. Gen zel et al.l 
2006), or that elliptical formation at high redshift might pro- 
ceed in more chaotic multiple mergers for which our simple 
approximations are not valid, especially in the most massive 
systems at very high redshifts z ~ 6. However, we have seen 
in Figure QT| that evolution in properties such as virial ve- 
locities, disk structure, and orbital parameters do not drive 
much additional evolutio n in M^/M^. Furthermore, some 
reassurance comes from iLi et al.l (120061) . who consider sim- 
ulations which adopt cosmologically-derived merger histo- 
ries for the brightest z ~ 6 quasars, involving multiple ma- 
jor and minor mergers, and find that the remnant obeys a 
similar M B h/M* relation to our idealized high-/ gas simula- 
tions - i.e. it is consistent with our estimated evolution in 
Mbh/M* (and R e /M*) as a function of redshift. Despite a 
chaotic, rapid merger history, the remnant properties can ap- 
proximately be predicted as a function of stellar mass and gas 
fraction from our simple prescription. The most important (at 
lowest order) element of our estimate is simply the qualita- 
tive statement that the progenitors of high-redshift ellipticals 
should be characteristically more gas-rich - which given the 
rapid star-formation timescales estimated for massive disks 
dBell & de J ong 2000) and observationally inferred brief, po- 
tentially burst-dominated star-formation hist ories of the most 
massive ellipticals Ce.g.. lThomas et afl l2005). is naturally ex- 
pected in theories of hierarchical growth. In other words, so 
long as the total stellar mass being violently scattered and to- 
tal gas supply reaching the spheroid center (i.e. total mass and 
gas fraction) are similar, the details of the merger history and 
distinction between binary or multiple mergers should not sig- 
nificantly change the simple gravitational physics which de- 
termine the most basic elements of the remnant structure. We 
therefore suggest that evolution in / gas is indeed the dominant 
(although not the only) physical agent driving evolution in the 
Mbh-M* relation. 

9. DISCUSSION 

Using a large set of numerical simulations of major galaxy- 
galaxy mergers, which include the effects of gas dissipation, 
cooling, star formation, and black hole accretion, we find that 
a feedback-driven model of BH growth and self-regulation 
predicts the existence of a BH "fundamental plane" (BHFP), 
of the form M BH oc a 30 R° e 5 or M BH ex M°/-° J a 15 ' 20 , anal- 
ogous to the FP of spheroids. Comparing with existing BH 
mass measurements, the observed systems appear to follow 
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a nearly identical BHFP relation. Specifically, there are sig- 
nificant (at > 99.9% confidence) trends in the residuals of the 
M^H-cr relation with M* and R e at fixed a, and likewise in the 
Mbh-M» relation (with a or R e at fixed M*). While changes 
in halo circular velocity, merger orbital parameters, progen- 
itor disk redshifts and gas fractions, ISM gas pressurization, 
and other parameters can drive changes in e.g. a at fixed M», 
and therefore change in the Mbh-ct or Mbh-M* relations, 
the BHFP is preserved. 

This provides a new paradigm for understanding the tra- 
ditional relations between BH mass and either bulge veloc- 
ity dispersion or mass. These correlations (as well as those 
with other bulge properties such as effective radius, central 
potential, dynamical mass, concentration, Sersic index, and 
bulge binding energy) are all projections of the same funda- 
mental plane relation. Just as the Faber-Jackson relation be- 
tween e.g. stellar mass or luminosity and velocity dispersion 
(M» - a) is understood as a projection of the more fundamen- 
tal relation between M*, cr, and R e , so too is the Mbh - o re- 
lation (Mbh oc cr 4 ) a projection of the more fundamental rela- 
tion Mbh oc a 3 R® 5 . Recognizing this resolves the nature of 
several apparent outliers in the Mbh _ cr relation, which sim- 
ply have unusual a values for their stellar masses or effective 
radii, and eliminates the strong correlations between residu- 
als (in both observations and simulations). While the various 
changes above in merger properties can and do bias the vari- 
ous projections of the BHFP to different values, they simply 
move remnants along the BHFP relation. 

Given the empirical tendency towards more compact, 
smaller-^ spheroids at fixed stellar mass M* at high redshift, 
the BHFP predict s that B Hs should be more massive at fixed 
M* . iTruiillo et alj ((2006) compile a number of observations of 
the sizes of early-type galaxies at fixed stellar mass (for typ- 
ical ~ galaxies), and find a best-fit trend R e oc (1 +z) . 
The observed BHFP predicts that BH mass scales roughly 
as oc M' 5 R~ l °, which yields the prediction that the typical 
hosted BH mass at fixed stellar mass (or ratio of Mbh/M*) 
should increase as (1 +z) 5 . This agrees well with recent di- 
rect e stimates of the BH to host stellar mass ratio at high red- 
shift dPeng et al.ll2006l) . as well as indirect estimates of the 
evolution in the mean Mbh/M* from comparisons of quasar 
lumino sity functions and ea rly-type mass density measure- 
ments dMerloni et al] 120041) . BH and spheroid mass func- 
tions dHopkins et al.ll 2006d). a nd quasar clustering as a func- 
tion of redshift dAdelberger & Steidel|l2005b IWvithe & Loebl 
2005; Hop kins et al.ll2007etlFine et alJl2006l) . Interestingly, if 
we consider this in greater detail, observations suggest that 
the evolution in sphe roid sizes is relatively weak to z ~ 0.8 
(iMcIntosh et alJ2005l) and stronger from z ~ 1 - 2. Our BHFP 
analysis argues that the same sho uld be true for the ratio of BH 
mass to stellar mass, and indeed IPeng et alj d2006l) note that 
there is no significant evolution at lower redshifts z < 1 in their 
sample, compared to the substantial evolution they observe at 
Z~l-3. 

We have also developed a physically motivated model for 
this evolution. Based on the empirical and theoretical expec- 
tation that the progenitor disks in typical mergers should be 
more gas-rich at higher redshifts, we expect mergers to be 
more dissipational, yielding more concentrated remnants and 
driving the evolution in Mbh/M* along the BHFP. Indeed, 
adopting an empirical estimate for the mean / gas as a func- 
tion of stellar mass and redshift, we predict a trend with red- 
shift in the size-mass relation of our merger remnants which is 



similar to the observations compiled in Trui illo et alj d2006l) . 
and consequently a trend in Mbh/M* like that observed by 
Peng et al. (2006). Our simulations thus provide critical sup- 
port to a rguments from a semi-an alytic context, such as those 
made by lKhochfar & Silk (2006), that the observed evolution 
in R e (M*) can be explained by the increasingly gas-rich, dis- 
sipational nature of merger progenitors at high redshifts. It 
is worth noting that, although it does not rule out such merg- 
ers occurring, the trend in R e (M*) can be explained entirely 
by changing gas fractions in gas-rich, dissipational mergers, 
without invoking subsequent "dry" mergers at low redshifts to 
increase R e (see also, Krause et al, 2007, in prepartion). 

We also emphasize that our results are en t irely consistent 
with the previous study of Robertson et al. (2006c). How- 
ever, in that case, the authors considered only the effects of 
the (relatively weak) scaling of disk sizes at fixed M. t with 
redshift, and found that this introduced a weak evolution in 
the Mbh _ cr relation. Allowing for / gas to scale systemati- 
cally with redshift drives the evolution in Mbh/M* which is 
analyzed herein, and placing both simulations and observa- 
tions in the context of the FP relation reco nciles the appar- 
ent disa greement between the predictions of Rob ertson et alj 
( 2006 c]) for the Mbh ~ cr relation and observations by e.g. 
IPeng et alj d2006l) of the high-redshiftMBH-M* relation. 

There are a number of direct, testable predictions of this 
fundamental plane model for the correlations between BH 
and host properties. At both low and high redshifts, systems 
should lie on the same BHFP. Therefore, m easurement s of the 
effective radii or velocity dispersions of the Peng et al. (2006) 
objects should find that they are more compact (smaller R e , 
larger cr) than their z = counterparts of the same stellar mass, 
in a manner consistent with the BHFP in Table Q] If it is re- 
ally the BHFP driving the apparent evolution in Mbh/M* with 
redshift, i.e. the fact that at fixed M*, higher-redshift systems 
are more compact, then this also predicts different evolution 
for BH mass relative to a (the Mb h ~ g relation) t han f or BH 
mass relative to M*. Adopting the iTrujillo et alJ d2006l) esti- 
mate for how R e scal es with redshift and the near-IR spheroid 
fundamental plane of lPahre et alj d 19981) to relate a and R e at 
fixed M*, along with our BHFP in terms of M* and cr, this 
predicts a trend of the form Mbh/c 4 oc (1 +z)~°' 25 , i.e. weaker 
and inverse evolution in Mbh/c at fixed stellar ma s s, quite 
similar to the predictions made by Robertson et Id] d2006ch 
and consistent with the observations of Shields et al. (2003). 

At low redshifts, improved measurements of the host prop- 
erties of systems with well-measured BHs can significantly 
improve constraints on the BHFP. As noted in Table [T] the 
present observations demand a correlation of the form Mbh oc 
a a M* over a simple correlation with either a or M* at > 3 cr 
confidence. Already, this puts strong constraints on theoreti- 
cal models of BH growth and evolution - BH mass does not 
simply scale with the star formation (stellar mass) or virial ve- 
locity of the host galaxy. However, there is still a substantial 
degeneracy between the slopes a and (3 (roughly along the 
axis (3 i=a 1 - a/4). For example, the current data do not allow 
us to significantly distinguish a pure correlation with spheroid 
binding energy Mbh oc (M* cr 2 ) 2 / 3 from the marginally fa- 

1/2 9 

vored relation oc M* a - both suggest that the ability of BHs 
to self-regulate their growth must be sensitive to the potential 
well at the center of the galaxy (and therefore to galactic struc- 
ture), but the difference could reveal variations in the means 
by which BH feedback couples to the gas on these scales. 
Increasing the observed sample sizes and, in particular, ex- 
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tending the observed baselines in mass and a will substan- 
tially improve the lever arm on these correlations. In partic- 
ular, the addition of stellar mass M* information to the sig- 
nificant number of objects which have measurements of a 
and indirect measurements of Mbh from reverberation map- 
ping would enable considerably stronger tests of our pro- 
posed BHFP relation. Furthermore, to the extent that the 
evolution in the Mbh/M* and R e (M*) relations is driven by 
the relatively gas-rich nature of the merger progenitors, the 
residuals in e.g. Mbh/M* should also be correlated with 
other tracers of the amount of dissipation in the spheroid- 
forming merger. These include quantities such as R e and er, 
of course, in the context we have discussed, but also e.g. 
the central phase-space density, kinematic properties such 
as rotation and kinematically decoupled c omponents (see 
iHernquist & Barneslll991t ICox et aLll200 6b). and potentially 
the presen ce of central cusps in the s tellar light profiles of the 
remnants (Mih os & Hernquis3ll994ah . We do note the caveat 
from § [7] however, that care should still be taken to consider 
only bulge properties and remove e.g. rotationally supported 
contributions to the velocity dispersion. 

Given the robust nature of the BHFP, we might also ask if 
there are processes which we might expect to drive systems 
away from the BHFP. For example, what are the effects of 
subsequent gas-poor (spheroid-spheroid or "dry") mergers on 
the BHFP? Such mergers, by definition, conserve total stellar 
mass and BH mass (simply adding the M» and Mbh of the two 
merged systems). However, simple energeti c arguments im- 
ply that a is not dramat i cally changed (e.g., Hernqui st et al.l 
I1993L iNipoti et al.ll2003b rBoylan-Kolchin et al.ll2006l) . If we 
assume a is unchanged in a "dry" merger, then the BHFP re- 
lations in Table[T]imply a ~ 0.08-0. 12 dex offset (in the sense 
of Mbh being too large) from a major (equal-mass) dry merger 



(the offset being ~ 0.03 -0.05 dex for a more probable 1:3 
mass-ratio merger). This is also supported b y a small subset 
of full numerical re-merger simulations from Robertson et al. 
(2006b). Although this implies that the BHFP will not survive 
a large number (> 3) of successive dry mergers, the observa- 
tionally estimated rate of ~ 0.5- 1 major dry mergers for a 
typical massive elliptical dBell et al.l2006Hvan Dokkuml2005l) 
implies that the realistic resulting deviations from the BHFP 
are small (smaller than the internal scatter in the relation it- 
self). This may be important, however, for explaining how 
the compact, high-redshift spheroids observed and predicted 
herein increase in size to lie on the local R e -M* relations 
(since each major dry merger will approximately double R e ), 
a possibility that will be investigated in detail in future work. 

Therefore, as appears to be borne out by the local obser- 
vations we consider, the BHFP appears to be a robust cor- 
relation, which provides an improved context in which to 
understand the nature and evolution of the numerous ob- 
served correlations between BH and host spheroid prop- 
erties. In particular, the results described here provide 
new, important constraints for models of BH growth, feed- 
back an d self-regulation, and support the proposal devel- 
oped by Hopki ns et al.l (l2006al) that major mergers between 
gas-rich galaxies represent the pr inciple mechanism for trig- 
gering intense starbursts ( e.g., Barnes & Hernquistl 1199 It 
Mihos & Hernquist 1994b, 1996) that evolve into quasars 
(e.g.. ISanders et all ll988) and which eventually leave rem- 
nants satisfying the sa me structural correlatio ns observed for 
elliptical galaxies (e.g., Robertson et al. 2006b). 
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ments. This work was supported in part by NSF grant 
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